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In this contribution we group the operator basis for d 2 dimensional Hilbert space in a way that 
enables us to relate bases of entangled states with single particle mutually unbiased state bases 
(MUB), each in dimensionality d. We utilize these sets of operators to show that an arbitrary density 
matrix for this d 2 dimensional Hilbert space system is analyzed by via d 2 + d + 1 measurements, 
d 2 — d of which involve those entangled states that we associate with MUB of the d-dimensional single 
particle constituents. The number d 2 + d+l lies in the middle of the number of measurements needed 
for bipartite state reconstruction with two-particle MUB (d 2 + 1) and those needed by single-particle 
MUB [(d 2 + l) 2 ]. 
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Two orthonormal vector bases, Bi, 
mutually unbiased bases (MUB) iff 

V \ui), \v,2) € B\, Bi resp. |(mi|m2)| = constant, 

(1.1) 

i.e., the absolute value of the scalar product of vectors 
from different bases is constant and independent of the 
vectorial labels within either basis. For a finite dimen- 
sion, d, Hilbert spaces the constant is Schwinger [1| 
was first to emphasize that there are more than two such 
bases "that exhibit maximum degree of incompatibility" 
i.e. more than just the pair of conjugate bases such as 
| a;), (spatial coordinates) and |fc) (momentum represen- 
tation basis). The information theoretic oriented term 
"mutually unbiased bases" (MUB) is due to Wootters Q. 
For the infinite dimensional Hilbert space case Wootters 
and coworkers 0, H, E| related the MUB's to lines in 
phase space. The transcription of these notions ("lines 
in phase space") to the finite, d, dimensional cases may 
be accomplished (0, 0, [H) via the eigenfunctions of the 
commuting operators 



x m Z^ x m z^ 



, I = si' 



with m,m',/,/',s integers. Here X and Z arc the 
Schwinger operators (SO) These are discussed in the 
following section, and abide by the relation ZX = uXZ, 
with u> = e % ~ . 

The maximal number of MUB was shown Q to be 
d+ 1. Ivanovic, [|| demonstrated that for d — p (a prime) 
the set of d + 1 MUB allows what is probably the most 
efficient means for determining the density matrix of an 
arbitrary state. The MUB analysis attracted great deal 
of research and cogent reviews are given, e.g., @, B @- 
These studies now involve abstract algebra and projective 
geometry: @, 0, EE EI El ES EE El d ■ Of particular 
interest to us are the articles by Planat and coworkers 
EH, EH who studied entangled states in conjunction with 
MUB sets similar to those studied in t he pres ent paper 
Similar ideas may be found also in 
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In this work we aim to relate e?-dimensional Hilbert 
space MUB (with d a prime) with bases of entangled 
states in d 2 -dimensional Hilbert space. To form this 
relation, we classify the basis for all operators in d 2 - 
dimensional Hilbert space via SO for d dimensional 
Hilbert space. Thus, we give a physical interpretation for 
the G? 2 -dimensional Hilbert space, as two particles each in 
d-dimensional Hilbert space. The basis for all operators 
in G? 2 -dimcnsional Hilbert space contains d 4 orthogonal 
operators. Following the rational of [f|, we cluster these 
operators into subsets of d 2 — 1 commuting operators, 
and construct their common eigenvalues. We show that 
these states form a basis in the c? 2 -dimensional Hilbert 
space. Each basis (correspond to different subset of op- 
erators) has a simple and transparent relation to onc- 
particle MUB. Some of the bases constitute maximally 
entangled states for the two-particle system. The re- 
lation between these entangled states and one-particle 
MUB is obtained via a projection of the one particle 
state onto the two particle entangled state. Finally, we 
apply our classification for the basis for all operators in 
<i 2 -dimensional Hilbert space for state reconstruction (to- 
mography). We show that one needs less measurements 
than are needed for bipartite state tomography analyzed 
only by single-particle MUB. 

In section II, utilizing formalism that stresses the (al- 
gebraic) field aspect 0] applied to the approach of 0, 
we obtain d + 1 MUB for the d dimensional one parti- 
cle system, that form subsets each of d — 1 commuting 
orthogonal operators. These operators are selected from 
the d 2 orthogonal complete set of operators that span 
the operator space. In this section cases with dimen- 
sionality, d, an odd prime are discussed. Extension to 
d powers of odd prime is discussed in appendix A. The 
analysis in section II aims to setup the stage for section 
III. In section III, we analyze the operator space of two- 
particle system, each belong to a <i-dimensional Hilbert 
space. In this section, we present the main results of this 
paper. First, we associate the single-particle d+ 1 sets of 
operators (each set contains d orthogonal commuting op- 
erators) with two-particle bases made of entangled states. 



2 



Then we show that the <i 2 -dimcnsional Hilbcrt space of 
two-particle density matrix is analyzed by d + d + 1 
measurements, d{d — 1) of which involve those entangled 
states that may be associated with single-particle MUB. 
We show that particularly symmetric density matrix may 
be accounted for by d 2 + 1 measurements, in close anal- 
ogy with the corresponding single particle MUB case. 
The numerology of the operator counting is given in Ap- 
pendix B. The last section includes conclusions and some 
remarks. 



II. MUTUALLY UNBIASED BASES 
DIMENSION 



FINITE 



Schwinger [l[ noted that the physics of finite dimen- 
sional, d, Hilbert space is expressed by via two unitary 
operators, X, and Z. Thus if we label the d distinct 
states, termed the computational basis, by |n), n = 
0, 1, • • • , d — 1; \n + d) = \n), these operators are defined 
by: 



Z\n) 



X\n) 



I) 



(2.1) 



with u> = e . They form a complete set, i.e. only a 
multiple of the identity commutes with both X, Z. We 
shall briefly outline a method to utilize these operators 
(due mainly to @) to construct the d + 1 MUB for a 
d dimensional Hilbert space with d being an odd prime. 
This review will be of help in building our sets of entan- 
gled states that we shall associate with these MUB. The 
computational basis vectors span the Hilbert space. All 
operators in this space are expressible in terms of the d 2 
Schwinger-operators [l[: 



X m Z l ; m,l = 0,1, 



.d- 1. 



(2.2) 



The operators X m Z l whose number is obviously d 2 form 
an orthogonal basis for all operators in the d dimensional 
Hilbert space, 



Tr 



X rn z l X m g 



— dS mm '5ll'. 



(2.3) 



This follows from Eq. (|2.ip which implies the commuta- 
tion formula 



XZ = uZX 



(2.4) 



Now, let us confine ourselves to cases wherein to, I £ F ( ; 
where is a Galois field with d elements. In this case, we 
can relate Schwinger operators to MUB. With this aim 
we group Schwinger operators (|2.2|) into d+ 1 sets of d— 1 
orthogonal commuting operators (which together with 
the identity operator form a complete operator basis for 
the Hilbert space) . Each set of (commuting, orthogonal) 
operators defines a unique vector basis in Hilbert space. 
All the d + 1 sets of bases form an MUB set. 

Let us first consider the case m = in Eq. (|2.2j) . Read- 
ily, the operators Z l with I = 0, 1, • • ■ , d— 1 (I = is, triv- 
ially, the identity operator) form one set of commuting 



and orthogonal operators. This set is diagonalized in the 
computational basis (c.f. Eq. (|2.ip ). Next consider the 
case m^O. In this case a unique inverse m _1 is defined 
on Fd and thus we can rewrite the operators (|2.2[) as 



X m Z l 



'(XZ") 



b\rn 



(2.5) 



where b = l/m, v = — %m(yci — 1) and b = 0, 1 ■ • • , d — 
1; to = 1, ••• ,d — 1. Thus, we have associated the 
d(d — 1) of the above operators in the following manner 



X m Z 



bur 



(xz b ) r 



(2.6) 



That is, these operators differ at most by a unimodular 
number. Now, for a fixed &, we have (d — 1) orthogonal 
and commuting operators. There are d distinct such sets, 
each labelled by b, b = 0, 1, • • • , d— 1 which arc orthogonal 
with 



Tr 



{xz b y 



(xz b 'y 



= d Sh.b'S. 



b,b' Qm.Tn' 



,to' ^ 0. 
(2.7) 

When each ^-labelled set of the d— 1 orthogonal, commut- 
ing operators is supplemented with the identity operator 
it constitutes a set of d unitary, orthogonal and com- 
muting operators. And as such, defines a vector basis 
for Hilbcrt space: The d (orthonormal) vectors diagonal- 
ize these operators in the set. (We remark, in passing, 
that the basis is defined up to a choice of a phase factor 
which does not affect the following results.) For each (b- 
labelled) set there exist a unique vector basis. Here, we 
are designating the vectors that form the basis by |6;c), 
where the index b labels the basis and c that particular 
vector in the basis b, (6, c = 0, 1, • • • , d — 1). The expres- 
sion for these states in terms of the computational basis 
is!, 



1 



E 

n=0 



n(n—l) — cn | 



(2.8) 



The eigenvalues are to c . 

The importance of the classification of the operators in 
this manner arises when considering the relation between 
different sets of the vector basis. One can readily check 
that these bases form MUB, i.e., 



(b;c\b]c') = S CtC >, 
\(b',c'\b,c)\ = b ^b>, 



(2.9) 



Thus, these d distinct sets of bases plus the computa- 
tional basis (which is mutually unbiased to all of these 
sets) form the maximal number of MUB, that is d + 1. 
The proof of the last formula involves the well known (l9j 
Gaussian sums. 

Since Schwinger operators, Eq. (|2.2p form an operator 
basis for d-dimcnsional Hilbert space, so do the set of 
operators (XZ b ) m together with Z b (b = 0, 1, ■ • ■ , d - 
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1; to = 1, 2, ■ • ■ , d— 1). Hence we may write an arbitrary 
density operator as 



/ d-l d-l 

-A EE Tr K^)1 ((xz b r) 



, m=l b=0 



d-l 



5>M (z l ) 



(2.10) 



We are now aiming to relate this form to the expression 
of density operators in terms of MUB's states \b;c). We 
shall first prove the following theorem 

Completeness Theorem 

Let several sets of operators, each labelled by b, con- 
stitute of d orthogonal, commuting unitary matrices, 
U a {b) (a = 0, 1, 2, • • • , d— 1), over a d dimensional Hilbcrt 
space. We take, without loss of generality, Uo(b) = I. Let 
\v b ), iy = 0, 1, 2, • • • , d — 1), be a basis that diagonalize 
the <i matrices of the & set, such that 



{v b \U a (b)\v' b ) = W?(b)5 v y, 
\W?(b)\ = 1. 

Then, the following (completeness) relation 

i£TO*)(W2(*))* = W 

holds. 
Proof 

Lets define d vectors for each set b, 



(2.11) 



(2.12) 



\Vb\Otb) 



1 



(2.13) 



Now the orthogonality of /7 a (6) implies 
dS a , a , = Tr [[/»(&)£&(&) 

= ^^(6)(^'(6))* = d^(a' 6 |^)(^|a fe ). 

:. (a'|a)=<W*'- (2.14) 

Thus the sets of d vectors \a b ) form a complete orthonor- 
mal basis for each b. Hence we have the completeness 
relation 

(u b \v' b ) = 8 Vb y b =J2(Mab)(a b H) (2.15) 

a 

□ 



Returning to the first term of Eq. (|2TTD|) - let (XZ b ) m 
correspond to U a (b), i.e. the index to corresponds to a 
(a > 0) and Uo(b) = I. Since the sum in first term of Eq. 
(|2.10p runs over to > 0, for each b, we add and subtract 
the identity. In the second term we have that Z l goes 
over to Z a - d orthonormal, commuting unitary matrices 
with Zq = I. With these replacements we may write, 

- d-l d-l 

b=o m=0 

-l 



rn— 



(2.16) 



Expressing each of these operators via its spectral repre- 
sentation, e.g., 



U m (b) 



d-l 

E 

c=0 



b;c)W™(b)(b;cl 



(2.17) 



where |6;c) are given in (|2.8p and substituting into the 
expressions above we get, after utilizing the completeness 
theorem, 



b.c 



\b;c)(b;c\p\b;c)(b;c\+Y / \^^\p\n)(n\ -I 

n 

(2.18) 



This is an expression for p in terms of probabilities. The 
number (6; c\p\b; c) corresponds to the probability to find 
the state \b;c) in p. We see that, as shown in Q, the 
state of the system, i.e. p, is obtainable via the d + 
1 measurements, where the measurement of a unitary 
operator is understood by the measurement of the two 
commuting hermitian operators, 



Mi = xz b + {xzy 
M 2 = i{XZ b - (xzy). 



(2.19) 



Each of these measurements yields the d — 1 independent 
probabilities (since the sum of the probabilities adds to 
one). This gives (d + l)(d — 1) = d 2 — 1 numbers that 
determine the density matrix. It should be noted that 
this holds since the operators are non-degenerate. 



III. ENTANGLEMENT ASSOCIATED WITH 
MUB 

We now wish to consider two d dimensional particle 
state. The Hilbert space dimensionality is now d 2 while 
the operator space is d 4 dimensional. Our aim is to uti- 
lize the single particle MUB given above to construct 
two particle entangled state bases that, thereby, may be 
viewed as associated with the MUB. 

We first wish to state what wc consider as a self evi- 
dent sufficient condition for entanglement: Given single 
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particle (particle n) operators A and A' which do not 
commute: AA' = otA' A, and similarly, given the non- 
commuting operators B and B' for a second particle (par- 
ticle v): BB' = (3B'B. Let a, (3 be scalars with a/3 = 1. 
The two two-particle operators AB and A'B' commute. 
We now assert that the common eigenfunction of these is 
an entangled state - we are unaware of a formal proof of 
this which we consider as obvious and use it as a guide 
for constructing entangled states. 

Returning to our two-particle Hilbert space. The dis- 
cussion of the previous section implies that the d A unitary 
orthogonal operators, 

(X mi Z l %(X m 'Z l ') v , m 1 ,h,m 2 ,l2 = 0A,--- ,d-l, 

(3.1) 

form a complete orthonormal operator basis. Restricting 
our consideration to mi,h, m2,h € IFd, allow us to re- 
construct and classify the d 2 operators in (|3.1[) as follows 



{XZ b % ll {XZ b2 )™ 2 , {d(d-l)) 2 operators 
z bi (xz b2 )™ 2 d-d(d-l) operators 
(XZ bl )™ 1 Z b v 2 d-d(d-l) operators 



Z b -Z b J 



d 2 operators , 



(3.2) 



here bi, b 2 = 0, 1, • ■ • , d — 1 and mi, 777.2 — 1, ■ • ■ , d — 1. 
This construction is carried out with an exact analogy to 
the one-particle case (c.f. Eq. (|2.5[1 ). One readily verifies 
that the above d 4 operators are orthogonal and thus span 
the operator space. 

First, let us note that these operators are diago- 
nalized in the one-particle MUB basis that was intro- 
duced in section II. Alas, these bases do not form an 
MUB set on the two-particle Hilbert space. Thus, for 
example, the basis which diagonalize Zj'flZ 1 ' 2 )"' 2 is 
{|n) M |&i; c)^}, while the operators Z b ^Z b v 2 are diagonal- 
ize in {\n') ^n) ,y} . These bases are clearly not MUB. 
However, in this construction we do find sets of MUB 
for the two-particle Hilbert space that is directly con- 
structed from the one-particle MUB. The basis which di- 
agonalize Z bl (XZ b2 )™ 2 (i.e., {\n)^\bi;c) v }) and the ba- 
sis that operators (AZ 61 )™ 1 ^ - {\bi; c) ^n) „} - are 
MUB. Another set of MUB are the basis that diagonal- 
ize {XZ b ^(XZ b2 )™ 2 ({|&i;c) M |6 2 ;c%}) and the basis 
that diagonalize Z b jZ b2 . 

We show that some of the above operators are diago- 
nalized by entangled states that form two-particle MUB 
and have a close relation to a single-particle MUB. 

To show this, we rewrite the first set in Eq. (|3.2[) as 



(x^z^uz^ry 1 , 

s, mi = 1, 2, • ■ • , d — 1, b, m = 0, 1, 2, • ■ • , d — (3.3) 

where the parameters s,m and b are defined as smi = 
777,2 ,sm = b\ and 62 = bm>2 — mm\. This set, for 
fixed s and b, contains d(d — 1) orthogonal, commut- 
ing operators. (Allowing s and b to vary, the sets contain 



[d(d— l)] 2 orthogonal operators.) Now consider d orthog- 
onal operators extracted from the last set of Eq. (|3.2p in 
the following way: We partition the set into two sets by, 
(r,« = l,2,-.-,d-l) 



z s ;z- r 



(d— l) 2 operators 
(d — 1) operators, 
(d — 1) operators, 



(3.4) 



these, with the identity operator, add up to the set of d 2 
orthogonal, commuting unitary operators. For fixed s, 
the (d— 1) operators Z^ r Z~ r commute with the d(d— 1) 
operators given in Eq. (|3.3[) defined above to form, with 
the identity, d 2 orthogonal commuting operators in our 
d 2 dimensional Hilbert. Thus these define the representa- 
tions wherein they are diagonal. In terms of the compu- 
tational basis for the two particles, these stats are given 
by 

\b,s; Cl ,c 2 )^ = -^Y, u,( ^ Mn ~ 1) ~ Cin \ n )»\ sn + C2 )»- 

* 71 

(3.5) 

Then it can be checked that (ci, c-x = 0, 1, 2 ■ ■ ■ , d— 1, we 
shall eschew the label fiis below) 

Z^Z~ 1 \b,s;c 1 ,c 2 ) = u~ C2 \b, s;c 1 ,c 2 }} 
{X) fl (X s Z sb ) v \b,s;c 1 ,c 2 ) = u c \b,s; Cl ,c 2 ), (3.6) 

c = sbc 2 . These states are entangled states (in accord 
with our guide), they are orthonormal and each set spans 
the d 2 dimensional Hilbert space. Since each operator 
set is labelled with the two indices, s and b, we have 
d(d — 1) distinct sets which constitutes a d 2 dimensional 
base composed of entangled states. Within these sets, 
d— 1 of them arc MUB. By direct calculation one obtains, 

\{b,s;ci,c 2 \b' ,s';c' 1 ,c' 2 )\ = 8 C1C2 b = b',s = s', 



Vd 



s', b^b', 



= _ s s > b b >. (3.7) 

Thus, bases with distinct s and b arc MUB (i.e. d-1 
entangled state bases) while common s and b defines a 
complete orthonormal d 2 dimensional basis. 

Next, we wish to relate these entangled state to a one- 
particle MUB. To convey the idea we focus on the set 
defined by s = 1, and without loss of generality we con- 
sider c 2 = 0: 



\b, s = 1; ci, c 2 = 0) M ,„ = -j= ^2 



u> 2 



n(n— 1)— c\n 



n),j\n) v . 
(3.8) 



Suppose that the /i particle is found in one of the states 
which belong to the MUB that were introduced in section 
II. Then, the reduced state for the v particle is (up to 
normalization) 

M (6i; c\b, s = 1; ci, c 2 = 0)^ = -^=|6-6i; a-c) v . (3.9) 
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This is the exact meaning of associating entangled states 
with single particle MUB. For the case b\ > b, the basis 
label for particle 2, b 2 = b — bi, is negative and should 
be defined as modulus d. Note that for an odd dimen- 
sion d, b\ is always different from b 2 . This ensures that, 
after measurement, states of the two particles belong to 
different MUB. The projection given in Eq. (|3.9p corre- 
sponds to Hcrmitian operators which were introduced in 
Eq. (I2~19l) . 

Finally, we gather the operators that provide a basis for 
the <i 2 -dimensional Hilbert space to state reconstruction. 
For this, we may use the d 4 unitary operators given in 
Eq. (f3~T]) as is rewritten via Eqs. (HH, (O, and (E3 to 
express an arbitrary two-particle density matrix as 



EE Tr 

b.s m,n 



x fl (x s z s \(z s u z- x y 



E [ tr p 



(tr P Z;(XZ m )t) (z2(XZ m )i\(3.W) 



E(^"K" 



I. (3.11) 



Using what we termed the "completeness" theorem that 
is proven in section II we rewrite this as (we eschew the 
particle number label viz. /1, v as being obvious in the 
expression below), where 

P = -1 + E E ( b is;c 1 ,c 2 \p\b,s;c 1 ,c 2 ) (3.12) 

b.S Cl,C2 

x \b,s;ci,c 2 ) fi ,v(b,s;ci,c 2 \ 
- (d- l)^(n,n'\p\n,n') x \n)^(n\ ® \n')„(n'\ 



+ EE^ ;c l^l^l")l 6;c ) x \n)^{n\®\b,c) v {b-c\ 

b=0 n,c 

+ ^^2(n\( b ;c\p\b;c)\n) x |6;c) M (6;c| <g> \n) v {n\ , 

b=0 n,c 

where s = 1, 2, • • • , d — 1, b,n,n' ,c\,c 2 = 0, 1, • • • , d — 1. 
In the above expression, p is given in terms of proba- 
bilities. Thus Eq. (|3.10p presents an operational scheme 
for two-particle state reconstruction which is based on 
Schwingcr operators of one-particle. For example, the 
probability of finding the system whose state is p in the 
(entangled) state \b,s;ci,c 2 ) two readings are required: 
that pertaining to measuring Z sb ) v and one that 

gives the value of measuring Z^Z' 1 . We consider these 
two readings as one measurement as it refers to the same 
sample. With this in mind an arbitrary density matrix 
is given in terms of d 2 + d + 1 measurements: The first 



term implies d{d — 1) measurements, the second one and 
the third and fourth d measurements each. The total 
number of measurements exceeds the optimal number, 
i.e., d + 1, since not all the bases in this construction 
are MUB. However, the number of measurements in this 
scheme is less than one would have needed by considering 
only single-particles MUB, that is, (d+ l) 2 = d 2 + 2d+ 1 
measurements. 



IV. CONCLUSIONS AND REMARKS 

Mutual unbiased bases for one particle states in d di- 
mensions may be characterized by states that diagonalize 
a set of d commuting orthogonal unitary operators which, 
in turn, is one of d + 1 such sets. The central attribute 
of these is that a measurement ascertaining the particle 
to be in a state \u) of one such basis implies that it to 
be with equal likelihood in any of the states \v) of any of 
the other bases. Maximally entangled state of two parti- 
cles central attribute is that partial tracing of the coordi- 
nates of one of the particles (i.e. performing non selective 
measurements on it) leaves the other particle with equal 
likelihood in any state. We relate these two notions by 
constructing d(d — 1) entangled states base vectors that 
form an MUB as well as relate to the one- particle MUB. 
We also gave an operational content for this interpreta- 
tion. We showed that an arbitrary two particle density 
matrix accounting for Hilbert space dimensionality d 2 , is 
accountable via d 2 + d + 1 measurements where d(d — 1) 
of the measurements involve entangled states which we 
related to the single particle MUB. 
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Appendix A: Operator count at finite dimensionality 

The numerology involves two counts. The first shows 
that our groupings of the operators retains all of them 
and each just once. The second involves the algebraic 
addition and subtractions of operators to obtain group- 
ings each of which involves d 2 — 1 (not including the iden- 
tity) orthogonal and commuting operators eligible for the 
application of the completion theorem, section II. We 
begin with the first demonstration: Eq. (|3.3[) contains 
d 2 (d — l) 2 orthonormal operators. Next we add the two 
sets of Eq. (|3.4[) that contain d 2 — 1 orthonormal opera- 
tors. Now we add the two sets that constitute the third 
and fourth lines of Eq. (|3.2p these these contain 2d 2 (d— 1 ) 
orthonormal operators. The sum total of these with the 
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identity added is d A . QED. 

The groupings of sets of d 2 — 1 of commuting (and orthog- 
onal) sets amendable for the applicability of the comple- 
tion theorem is done as follows: first we add to each of the 
d(d— 1) commuting operators labelled by distinct b and s 
the set of d — 1 operators (Z^ l Z~ 1 ) n , n = 1, 2, • • • , d — 1. 
The combined set is composed of d 2 — 1 orthogonal 
and commuting operators (to which, when supplemented 
with the identity, we associate the d 2 dimensional or- 
thonormal basis \b, s; ci, C2)). To compensate we sub- 
tract this added set and then combine this subtracted 
set with the rest of the operators. Now we add to it 
— (d— 1) J2 n =i {ZJl + Z™) ' these together (supplemented 
with the identity) form d — 1 sets each with d 2 commut- 
ing and orthogonal unitary operators. Correcting for this 
subtraction we add to each of d sets {b = 0, 1, • • • , d— 1) 
of the form J2 n =o m=i ^(^^ b )"' to § ct a g am sets each 
of which is made up (with the identity) of d 2 commuting 
and orthogonal unitary operators. 



Here b,c,n 6 F<j, is Galois field with d elements, 
|n) vector in the computational basis (labelled with an 
element of the field); and tr[a] = a+a p +a p +■ ■ -+a p 
The trace ,tr, is a mapping with a g F^; trae¥ p . Now 
a basic property of trace is: tr[a+/3] = tr[at]+tr\]3]. (The 
factor i in the exponent means the solution of 2x = 1 6 
F,j.) Using our entanglement scheme (cf. fl5l . (HI), i.e. 
entangling within the computational basis, the entangled 
state corresponding to the generic MUB state Eq. (14. ip 
is, 



|6;c)^ = -^5>f [| " 2+m] |n>>>„ (4.2) 

neF d 

1 trf^r Ti 2 +cinl tr[^3-n 2 +C2nl I . i > 

neF d 



Appendix B: Entanglement of MUB states at finite 
dimensionality 

Our method for the finite dimensional may be applied 
to the case where the states' label are finite field variables. 
This is the case for dimensionality d, with d = p" , p an 
odd prime p], OJ, EE Ell ■ The MUB are S iven b Y> usin S 
our notation, (cj p = e l p ): 



\b;c) 



Vd ^ 



tr[^n 2 +cn] . 



(4.1) 



with bi + 62 = ci + C2 = c. Thus projecting the 

first particle to (any of the) 61, C\ labelled state results 
in projecting the second particle to \b 2 , c-i) v state. 



( (6i;ci|6;c) MI/ = -j=\b 2 ; c 2 ) v . (4.3) 



new d 
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